Spatial solitons in nonlinear left-handed metamaterials 
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We predict that nonlinear left-handed metamaterials can support both TE- and TM-polarized self- 
trapped localized beams, spatial electromagnetic solitons. Such solitons appear as single- and multi- 
hump beams, being either symmetric or antisymmetric, and they can exist due to the hysteresis-type 
magnetic nonlinearity and the effective domains of negative magnetic permeability. 
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Left-handed metamaterials, i.e. materials with simul- 
taneously negative real parts of dielectric permittivity 
and magnetic permeability |1J, demonstrate many pe- 
culiar properties including negative refraction and sub- 
wavelength imaging. A number of recent experiments 
reported the observation of some of those properties 
for artificially fabricated composites 0, 0] and planar 
transmission-line structures Q at microwave frequencies. 
There are strong efforts to fabricate the left-handed com- 
posite materials operating at other frequencies, and neg- 
ative magnetic permeability at THz frequencies has been 
recently demonstrated in planar structures composed of 
nonmagnetic conductive resonant elements p|- 

Composite materials used for demonstrating the prop- 
erties of left-handed media and negative refraction pos- 
sess combined negative electric and magnetic response in 
some finite frequency range determined by the resonant 
properties of nonmagnetic conductive elements of the ar- 
tificial composites. Nonlinearity introduced in such ma- 
terials would allow to control the frequency domain with 
the left-handed properties through strong modification 
of the effective nonlinear magnetic response 0. In par- 
ticular, both dielectric permittivity and magnetic perme- 
ability of the composite material consisting of a mesh of 
wires and split-ring resonators embedded into a nonlin- 
ear dielectric can be controlled by varying the electro- 
magnetic field intensity in such a way that it may allow 
switching the material properties between left- and right- 
handed [3 • Since microscopic fields can be enhanced dra- 
matically in a composite structure due to the split-ring 
resonators, strong nonlinear effects arc expected to be im- 
portant in such types of metamaterials, and this can stim- 
ulate applications of nonlinear left-handed materials for 
switching devices, limiters, and frequency conversion. In 
addition, different types of nonlinear elements can also be 
easily introduced into planar left-handed transmission- 
line structures Q making a response of such a structure 
strongly nonlinear. 

In this Letter, we study, for the first time to our 
knowledge, self-trapping of electromagnetic waves and 
the formation of spatial solitons in nonlinear left-handed 
composite media, and find the conditions for the exis- 
tence of both TE- and TM-polarized bright and dark 
solitons @- We demonstrate that left-handed materials 
with a hysteresis-type (multi-stable) nonlinear magnetic 



response support novel and unique types of single- and 
multi-hump (symmetric, antisymmetric, or even asym- 
metric) backward-wave spatial electromagnetic solitons 
due to the effective domains corresponding to different 
branches of nonlinear magnetic permeability. 

We consider a nonlinear composite structure described, 
within an effective medium approximation, by the non- 
linear magnetic permeability of the form [6j 
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where H is the applied magnetic field, F is a filling factor 
of the composite, ojqnl is a nonlinear eigenfrequency of 
the split-ring resonator, and V is the loss coefficient. We 
assume that the infilling dielectric of the structure has a 
Kerr-like nonlinear response, 



e D (\E\ 2 ) = e D0 + a\E\ 2 /E 2 c , 



(2) 



where cdo is the linear dielectric permittivity, E c is a 
characteristic electric field of a nonlinear response, and 
a = ±1 stands for focussing or defocusing nonlinearity, 
respectively. Here we neglect the dependence of the effec- 
tive dielectric permittivity on the electric field assuming 
that a nonlinear dielectric fills only the slits of the split- 
ring resonators, as discussed in Ref. [(|. 

Dependence of the nonlinear eigenfrequency of split- 
ring resonators on the macroscopic magnetic field H can 
be found in an implicit form (see Eq. (11) of Ref. |||), 
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where X = ujqnl/luq is the normalized nonlinear fre- 
A 2 = 16e 3 D0 LO 2 h 2 / c 2 , Q = uj/uq, uj = 



quency, 



(cj a)[d g /irhtDo] 1 / 2 is the linear frequency of a split-ring 
resonator, L' = T/loq is the normalized damping coeffi- 
cient, c is the speed of light, d g , h, a are the resonator 
parameters, and the magnetic field H is normalized to 
the critical electric field E c . 

When the losses are negligible, the dependence of the 
magnetic permeability on the normalized magnetic field 
H can have several distinct forms [6j . The most interest- 
ing, multi- valued dependence is shown in the inset Fig. ^ 
for the defocusing nonlinearity and Q > 1. 
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FIG. 1: (color online) Phase trajectories on the plane 
(H, dH/dx) described by Eq. Q at 7 = 2.8 with the nonlinear 
magnetic response shown in the inset for a — — 1, too = 12.8, 



h 



0.03 cm, d g 



0.01 cm, a 



0.3 cm, and F 



0.4. 



Solid, dashed, and dotted lines mark the phase trajectories 
corresponding to the solid, dashed and dotted branches of 
the nonlinear magnetic permeability. 



TM-polarized solitons. First, we study spatially lo- 
calized TM-polarized waves that are described by one 
component of the magnetic field and two components of 
the electric field. We consider monochromatic station- 
ary waves with the magnetic field component H = H y 
propagating along the z-axis and homogeneous in the 
y-direction, [~ exp (iwt — ikz)], described by the dimen- 
sionless nonlinear Helmholtz equation 



d 2 H 

dx 2 



+ [e^\H\ 2 )- 1 2 }H = 0, 



(4) 



where 7 = kc/uj is a wavenumber, x = x'lo/c is the di- 
mensionless coordinate, and x' is the dimensional coor- 
dinate. Different types of localized solutions of Eq. @ 
can be analyzed on the phase plane (H, dH/dx) (see, 
e.g., Refs. First, we find the equilibrium points: 

the point (0, 0) existing for all parameters, and the point 
(0, Hi), where Hi is found from Eq. (J3J at 



x 2 (Hi) = x 2 =nUi 



Fe cS 
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Below the threshold, i.e. for 7 < 7 tr , where 7^. = 
e[l + FQ 2 /(1 - O 2 )], the only equilibrium state (0,0) 
is a saddle point and, therefore, no finite-amplitude or 
localized waves can exist. Above the threshold value, i.e. 
for 7 > 7t r , the phase plane has three equilibrium points, 
as shown in Fig. ^ for the particular value 7 = 2.8. 

In the vicinity of the equilibrium state (0,0), linear 
solutions of Eq. @ describe either exponentially grow- 
ing or exponentially decaying modes. The equilibrium 
state (0,-Hi) describes a finite-amplitude wave mode of 
the transverse electromagnetic field. In the region of 
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FIG. 2: (color online) Examples of different types of solitons: 
(a) fundamental soliton; (b,c) solitons with one domain of 
negative or positive magnetic permeability (shaded), respec- 
tively; (d) soliton with two different domains (shaded). Insets 
in (b,c) show the magnified regions of the steep change of the 
magnetic field. 



multi-stability, the type of the phase trajectories is de- 
fined by the corresponding branch of the multi-valued 
magnetic permeability. Correspondingly, different types 
of the spatial solitons appear when the phase trajecto- 
ries correspond to the different branches of the nonlinear 
magnetic permeability. 

The fundamental soliton is described by the separa- 
trix trajectory on the plane (H, dH/dx) that starts at 
the point (0,0), goes around the center point (0,Hi), 
and then returns back [see Fig. solid curve]; the cor- 
responding soliton profile is shown in Fig. EJa). More 
complex solitons are formed when the magnetic perme- 
ability becomes multi-valued and is described by several 
branches. Then, soliton solutions are obtained by switch- 
ing between the separatrix trajectories corresponding to 
different (upper and lower) branches of magnetic per- 
meability. Continuity of the tangential components of 
the electric and magnetic fields at the boundaries of the 
domains with different values of magnetic permeability 
implies that both H and dH/dx should be continuous. 
As a result, the transitions between different phase tra- 
jectories should be continuous. 

Figures|2b,c) show several examples of the more com- 
plex solitons corresponding to a single jump to the lower 
branch of n(H) (dotted) and to the upper branch of n(H) 
(dashed), respectively. The insets show the magnified 
domains of a steep change of the magnetic field. Both 
the magnetic field and its derivative, proportional to the 
tangential component of the electric field, are continu- 
ous. The shaded areas show the effective domains where 
the value of magnetic permeability changes. Figure Hfd) 
shows an example of more complicated multi-hump soli- 
ton which includes two domains of the effective magnetic 
permeability, one described by the lower branch, and the 
other one - by the upper branch. In a similar way, we 
can find more complicated solitons with different number 
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FIG. 3: (color online) Normalized soliton energy flow vs. 
its propagation constant. Vertical solid line is a boundary 
of wave localization in the linear limit, j = 7t r . Thin verti- 
cal dashed line shows the existence boundary for the funda- 
mental soliton. Thick solid curve shows the dispersion of a 
fundamental soliton. Dashed and dotted curves correspond 
to the solitons with one domain where the value of magnetic 
permeability changes. Dotted curves correspond to one-hump 
solitons shown in Fig. I2lbl , dashed curves correspond to the 
two- hump solitons shown in Fig. EJ C ) - 



of domains of the effective magnetic permeability. 

We note that some of the phase trajectories have dis- 
continuity of the derivative at H = caused by infinite 
values of the magnetic permeability at the correspond- 
ing branch of fi c ff(H). Such a non- physical effect is an 
artifact of the lossless model of a left-handed nonlinear 
composite considered here for the analysis of the soliton 
solutions, tn more realistic models that include losses, 
the region of multi-stability does not extend to the point 
H = 0, and in this limit the magnetic permeability re- 
mains a single-valued function of the magnetic field 

For such a multi- valued nonlinear magnetic response, 
the domains with different values of the magnetic per- 
meability "excited" by the spatial soliton can be viewed 
as effective induced left-handed waveguides which make 
possible the existence of single- and multi-hump soliton 
structures. Due to the existence of such domains, the 
solitons can be not only symmetric, but also antisym- 
metric and even asymmetric. Formally, the size of an ef- 
fective domain can be much smaller than the wavelength 
and, therefore, there exists an applicability limit for the 
obtained results to describe nonlinear waves in realistic 
composite structures. 

Nonlinear dispersion of several types of spatial solitons 
is shown in Fig.|31as the dependence of the soliton energy 
flow vs. the normalized propagation constant 7. The fun- 
damental soliton exists in a limited range of 7. The lower 
boundary is determined by the threshold value 7t r . At 
the upper boundary, the amplitude of the soliton field ex- 
ceeds the critical value H* , and the magnetic permeabil- 
ity becomes positive. All spatial solitons are backward 
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FIG. 4: (color online) Phase diagram of Eq. (|SJ corresponding 
to the multi- valued magnetic response shown in Fig. Pf inset). 
Solid, dashed and dotted curves mark the phase trajecto- 
ries that correspond to solid, dashed and dotted curves in 
Fig. 0^inset), respectively. Inside the circle with the radius 
H* the magnetic permeability is a three-valued function of 
H, but it is single- valued, otherwise. Thick solid figure-eight 
curve is defined by a singularity in © , at which the (dotted) 
phase trajectories change their direction. 



propagating waves here, since the energy propagates in 
the direction opposite to the wavevector, as indicated by 
the negative sign of power flow. 

When the infilling dielectric of the structure displays 
self-focusing nonlinear response, we have 51 < I, and in 
such system we can find dark solitons, i.e. localized dips 
on the finite-amplitude background wave Q. Similar to 
bright solitons, there exist both fundamental dark soli- 
tons and dark solitons with domains of different values 
of magnetic permeability. For self-defocusing nonlinear- 
ity and Q < I, magnetic permeability is a single- valued 
function, and such a nonlinear response can support dark 
solitons as well, whereas for self-focusing dielectric, and 
51 > I, no solitons can exist. Thus, in the composite with 
self-defocusing dielectric, depending on the frequency of 
the electromagnetic field, either bright or dark solitons 
can exist. 

TE-polarized solitons. The TE-polarized waves are de- 
scribed by one component of the electric field and two 
components of the magnetic field. For such waves, mag- 
netic permeability depends on two components of the 
magnetic field and the analysis of localized modes be- 
comes more involved. 

We look for the waves propagating in the z— direction, 



E = y Ee 



iujt—ihz 



H = {-x H x + z iH z y- 



where (xn,yo,Zo) are the unit vectors, and E, H x , and 
H z are scalar real amplitudes of the fields. A constant 
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FIG. 5: (color online) Energy flow of the fundamental soliton 
(solid) and two solitons with a single domain (dashed and dot- 
ted). Solitons characterized by the dashed curves (left inset) 
have a single domain (shaded) with magnetic permeability 
described by the dashed curve in Fig.^inset). Solitons char- 
acterized by the dotted curves (right inset) have a domain 
(shaded) with magnetic permeability described by the dotted 
curve in Fig. H]a, inset). Examples show the profiles of the 
components H x (solid) and H z (dashed) in the TE soliton. 

phase difference between the components E and H z is 
assumed, and no energy flow in the x-direction occurs. 
Spatial distribution of the magnetic field in a soliton is 
described by a system of coupled equations, 

^ = [ £M (H 2 )- 7 2 ]^/ 7 , 

dH* 7 2 /T(g 2 )+C[^(g 2 )-7 2 ] rfi , 
dx [7M(# 2 )+7£] ' [ ' 

where £ = 2H 2 f/(H 2 ), H 2 = H 2 + H 2 . We analyze 
Eqs. © on the plane (H X ,H Z ), as shown in Fig- El 

Similar to the TM-polarized waves, the phase trajec- 
tory on the plane (H X ,H Z ) is defined by a particular 
branch of the multi-valued nonlinear magnetic perme- 
ability. Near the origin (0,0), the phase trajectory is 



defined by the solid branch of the magnetic permeabil- 
ity [see Fig. [IJ inset)]. In general, the phase diagram has 
three equilibrium states, which are similar to those de- 
scribed above for the TM-polarized waves. A singularity 
curve (thick solid line) is defined by the vanishing de- 
nominator in the second equation of Eqs. ©. The phase 
trajectories represented by dotted curves change their di- 
rection when they cross the singularity line. 

For a homogeneous structure with the nonlinear mag- 
netic permeability shown by a solid curve in Fig.^inset), 
the separatrix curve on the phase plane of Fig. 0] de- 
scribes the fundamental two-component TE-polarized 
spatial solitons with a smooth envelope. However, there 
exist more complex TE-polarized solitons, and such soli- 
tons include effective domains with different values of the 
multi-valued magnetic permeability. At the boundaries 
of the domains, both the fields H x and H z //j, should be 
continuous. As a result, the phase trajectories describ- 
ing the TE-polarized solitons with such domains will be 
discontinuous on the plane (H X ,H Z ). In Fig. we show 
the energy flow for the family of the fundamental soli- 
tons and two types of solitons with a single domain, as 
the functions of the normalized propagation constant 7. 
Example of the solitons presented here correspond to a 
symmetric dependence of the x-component of the mag- 
netic field; however, both antisymmetric and asymmetric 
TE-polarized solitons can exist in the structure as well. 

In conclusion, we have analyzed both TE- and TM- 
polarized self-trapped nonlinear localized beams - spatial 
electromagnetic solitons - in left-handed metamaterials 
with nonlinear resonant response. We have revealed the 
existence of unique types of solitons supported by the 
hysteresis-type magnetic nonlinearity and the domains 
with different values of magnetic permeability. Such soli- 
tons can have symmetric, antisymmetric, or asymmetric 
structure. We believe that similar solitons can be found 
in other types of complex nonlinear materials with left- 
handed properties or the frequency-dependent domains 
corresponding to negative refraction. 

We thank Alex Zharov for useful discussions and ac- 
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